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We consider the superconducting vortex solid-to-liquid transition in heavy ion-irradiated un- 
twinned YBa2Cu307_i single crystals in the case where the magnetic field direction does not coin- 
cide with that of the irradiation-induced linear columnar defects. For a certain range of angles, the 
resistivities measured in three orthogonal spatial directions vanish at the transition as three clearly 
different powers of reduced temperature. At previously known second-order phase transitions, scal- 
ing of physical quantities has either been isotropic or anisotropic in one direction. Thus, our findings 
yield evidence for a new type of critical behavior with fully anisotropic critical exponents. 



Vortex matter in high temperature superconductors 
offers many opportunities to study the effects of disor- 
der, fluctuations, and frustration, which lead to a va- 
riety of new phases and phase transitions Among 
these, the two principal thermodynamic phases are a dis- 
sipative vortex liquid at high temperatures, and a truly 
superconducting vortex solid at low temperatures. In 
clean superconductors, the vortex solid-to-liquid transi- 
tion is a first-order melting transition. The presence of 
strong disorder is believed to turn it into a continuous 
isotropic "vortex glass" transition ■ Then, many quan- 
tities develop power-law singularities upon approaching 
the transition, which define its critical exponents. The 
critical properties are insensitive to microscopic details 
such as small-scale anisotropies, and depend instead on 
such general features as symmetry and dimensionality, 
which determine the transition's universality class. All 
phase transitions discovered so far are of two main types, 
those with isotropic scaling properties, and those with 
anisotropic scaling in one direction. 

For point disorder the glass transition will be isotropic, 
i.e., the critical exponents are independent of direction, 
even though the material itself can be anisotropic, as is 
the case in, e.g., YBa2Cu307_a (YB CO). A markedly dif- 
ferent situation occurs in superconductors containing lin- 
ear columnar defects following irradiation by swift heavy 
ions 0, 0, S @, 0, H H- Such defects are very effec- 
tive at pinning the vortex lines when a magnetic field is 
applied parallel to them; the vortex solid at low temper- 
ature is then a "Bose glass" Q. The columnar defects 
and the applied field both break the symmetry, making 
the Bose glass transition anisotropic, with different ex- 
ponents parallel and perpendicular to the columns 
An interesting situation occurs when the magnetic field 
is tilted away from the defects. For small tilt the vor- 
tices will stay on the columns, leading to a transverse 
Meissner effect. As the tilt angle or the temperature is 
increased the vortices will eventually depin and enter the 



vortex liquid phase. The transverse component of the 
magnetic field then makes all three directions nonequiv- 
alent and opens for the possibility of fully anisotropic 
scaling, i.e., with different critical exponents in all direc- 
tions [10| | . Furthermore, the problem of vortex depinning 
from columnar defects is formally equivalent to the Bosc- 
glass transition of bosons in a random potential @|. In 
this analogy, the tilted magnetic field corresponds to an 
imaginary vector potential, resulting in a non-Hcrmitian 
localization problem (ll| . 

Here, we investigate the vortex solid-to-liquid transi- 
tion by electrical transport measurements performed in 
three orthogonal spatial directions on untwinned heavy- 
ion irradiated YBCO with the magnetic field tilted away 
from the columns. It appears that, for angles not too 
close to the principal crystal axes, the resistivity vanishes 
as a power law of reduced temperature with different ex- 
ponents in the three directions. Thus, we provide exper- 
imental evidence for a phase transition with anisotropic 
critical exponents in all directions, thereby extending the 
possible type of universality classes. 

YBCO single crystals were grown by a self-flux method 
in yttria stabilized zirconia crucibles As-grown un- 
twinned crystals with typical sizes of 0.5 x 0.2 x 0.02 mm 3 
were annealed at 400 °C in flowing O2 for a week. Such 
crystals show clear vortex melting signatures in mag- 
netic fields. Irradiation with 1 GeV Pb 56+ ions was 
performed at the Grand Accelerateur National d'lons 
Lourds (GANIL), Caen, France, with the incident beam 
almost parallel to the crystallographic c-axis. The result- 
ing damage consists of randomly distributed amorphous 
tracks, with a diameter of about 7 nm, extending through 
the total thickness of the sample 13]. Every ion impact 
creates a track. For most samples, the track density was 
iid = 1 x 10 11 ions cm' 2 , corresponding to a matching 
field B,), = 4>o/nd = 2.0 T ( tpo = h/2e is the flux quan- 
tum). Electrical contacts were prepared by silver paint, 
giving contact resistances below 1.5 tt. For resistivity 
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FIG. 1: (Color online). Experimental geometry. The colum- 
nar defects are parallel to the z direction, which coincides with 
the crystallographic c axis. The magnetic field is always ap- 
plied in the rrz-plane, at an angle 6 with respect to z. Vortex 
lines may accommodate to the columns by forming kinks. 



measurements in the x and y directions, perpendicular 
to the columns and parallel to the Cu02 planes, the cur- 
rent pads were applied on opposite sides of the crystal. 
The resistivity parallel to the columns, p z , was measured 
on another sample, with current contacts covering most 
of the ab-planc surfaces so as to assure the most homo- 
geneous current distribution 14]. The linear resistivity 
was measured using the standard four-probe technique, 
with a current of 0.1 mA and 1 mA for the in-plane and 
c-axis measurements respectively. Using a dc picovolt- 
meter as preamplifier, the voltage resolution was below 
1 nV. The angular resolution of the field orientation is 
0.01° along the direction of rotation of our single-axis 
rotating sample holder and 1-2° perpendicular to it. 

The geometry of our experiment is shown in Fig. [1] 
The magnetic field, H , is always applied in the xz plane 
at an angle 9 with respect to the direction of the columnar 
defects; p x>y were measured on the same sample mounted 
in different orientations with respect to the field direc- 
tion. Fig. [2] shows raw p(T) data measured in the three 
spatial directions. The zero- field transition temperature 
for the crystal used for most of the p x and p y measure- 
ments was T cQ = 91.7 K (AT C w 0.5 K); the crystal 
employed to measure p z had T c o = 91.2 K (AT C w 1 K). 
The difference in T c o accounts for the different vortex 
solid-to-liquid transition temperature T c found in the p z 
measurements. Our results are independent of the choice 
of x as either the crystallographic a or & axis. 

Close to the vortex solid-to-liquid phase transition, the 
resistivity drops to zero as a power law of \T/T C — 1|, 
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From the measured p%(T) curves, Si and T c were deter- 
mined in the standard way. A straight line was fit to the 
inverse logarithmic derivative of the experimental data, 
{dXnp./dT)" 1 = (T - T c )/s h as shown in Fig. H Be- 
low 5% of the normal state resistivity, power-law scal- 
ing corresponding to the critical regime is observed. All 
measurements are made in the linear response regime; 
therefore, the scaling is insensitive to circumstances such 
as possible non-locality of vortex flow. We have also di- 
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FIG. 2: (Color online). Normalized resistivity in three in- 
dependent directions of an untwinned heavy-ion irradiated 
YBCO single crystal with a matching field = 2 T, and 
at tilt angles (a) 9 — 4° and (b) 6 = 45° from the columns. 
Dashed lines are fits to Eq. (JTJ. Inset: Extraction of T c and 
the critical exponent s; using (<91n pi/dT)' 1 = (T — T c )/si. 



rectly fitted the data to a power law (see Fig. \2§ , which 
yields consistent results. 

The experimentally obtained exponents s x , s v and s z 
for poH = B<f, are shown as function of 9 in Fig. [Ha). 
There are three distinct angular regimes, with different 
relations between the three exponents. For 9 < 15°, s x ~ 
s y , while s z is smaller. For 15° % 9 < 65° all Si are 
different, with s y < s x < s z . In this interval the best 
fits to the power law |T]) were found. The exponents are 
approximately constant as a function of angle, suggesting 
a certain amount of universality. The average values s x = 



4.9 ± 0.2, s v = 3.6 ± 0.2, and s z = 6.9 ± 0.3 are indicated 
by the dashed lines in Fig. Ela) . Finally, at angles above 
65°, s z decreases while s x increases. The error bars were 
estimated by varying the resistance interval of the fits 
while maintaining good agreement with the data. 

The interaction between the columns and the vortices 
changes with tilt angle, which can explain the existence 
of several angular regimes with different behavior of the 
exponents. At angles below 9l ~ 1 — 2° [lij], the vortices 
are locked within the columns in the Bose-glass phase and 
a transverse Meissner effect is expected 0, UH . The prox- 
imity to this phase will influence the transition and may 
explain the variation in the exponents for 9 < 5° [3]. For 
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the sequence of the exponents is observed, probably due 
to the layered nature of YBCO 
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FIG. 3: (Color online). Exponents s x , s y and s z of the power- 
law Eq.fTJ by which the resistivity vanishes at the vortex 
solid-to-liquid transition in the x, y and z directions respec- 
tively, (a) Exponents at 2 T as function of angle 9. For 
15° < i$65°, the scaling behavior is anisotropic in three di- 
rections, (b) Exponents at 8 = 45°, as function of magnetic 
field. There is anisotropic scaling at fields close to the match- 
ing field, noH sa = 2 T. Dotted lines show the average 
value of Si in the regimes of anisotropic scaling. The error 
bars are estimated from the quality of the fits in FigO 



9 > 6l, vortices adjust to the columnar defects by form- 
ing a staircase structure 0, Q , depicted in Fig. [Q This 
occurs up to a certain accommodation angle 9 a deter- 
mined by the balance between the pinning energy gained 
by the segments trapped on a column and the energy 
loss due to vortex line tilt. For 9 > 8 a , the columns are 
thought to act as translationally symmetric "point pins" . 
The vortex lines intersect the columns only over a length 
corresponding to the vortex core radius, and are on aver- 
age parallel to the direction of the applied field. Exper- 
imentally, 9 a is usually defined as the maximum of the 
angular-dependent resistivity I(| 17 1 , as shown in the in- 
set of Fig.0J In our experiments this gives 9 a « 15 — 20°, 
which coincides with the observed change in the behav- 
ior of the critical exponents. However, it is above this so 
defined 9 a that three different resistivity exponents are 
observed. This suggests that the columnar defects still 
play a determining role in this regime, as they break the 
symmetry sufficiently strongly for the anisotropic scaling 
to be observed. Finally, as the field orientation perpen- 
dicular to the columns is approached, another change in 



Fig. [3^b) shows that the values of Sj change as a func- 
tion of H , indicating a change in the nature of the phase 
transition. At high fields, the vortex density exceeds the 
defect density. Hence, vortex- vortex interactions become 
increasingly important compared to the vortex-defect in- 
teraction, and the anisotropy of the transition described 
above should progressively disappear. Above 4 T, the ex- 
ponents are roughly the same in all directions and more 
or less field independent. The difference between the ex- 
ponents is most pronounced close to the matching field. 
The average values in this region are s x — 4.7 ± 0.2, 
s y = 3.8 ± 0.2 and s z = 6.6 ± 0.3, consistent with the av- 
erage values at constant ^qH = 2 T found above. At 
lower fields, the exponents again approach a common 
value. We observe these trends systematically irrespec- 
tive of sample or matching field. 

For high fields the effect of the columns is less rele- 
vant. Most vortices are pinned by pointlike pinning cen- 
ters, or by the interaction with the "matrix" of vortices 
trapped on a column [lfi ] . The vortex solid resembles an 
isotropic glass, with identical power-law exponents in all 
directions. Then, the glass line is given by [2pj ]. 



H g (T) oc (1 - T/T c0 ) n , 



(2) 



where T c q is the zero-field transition temperature and 
H g = H(T = T C {H)). Fig. H shows the experimentally 
obtained solid-to-liquid transition temperatures T C (H) 
obtained from p y . To account for the intrinsic anisotropy 
of YBCO an effective field H cS = H^~r 2 sin 2 9 + cos 2 9 
is plotted on the ordinate ( 7 w 7 is the anisotropy pa- 
rameter 2l|). The Figure also shows a fit to Eq. ([2]) of 
the T c data at fields where approximately isotropic expo- 
nents St were found. The fit gives rt ~ f.5, in agreement 
with the value reported for the isotropic vortex glass [2(| • 
We observe that for fields close to the matching field, 
and angles 9 < 65°, T C (H) lies above the glass line (|2|). 
A bend towards this line is observed at an angle consis- 
tent with the experimentally defined 9 a ss 15° but, up to 
9 r~j 65°, the data still lie above the H g (T) line. In fact, 
T C (H) resembles the behavior expected for the Bose-glass 
line [f , 22|, 23| , adding further evidence that the effect of 
the columns as correlated pinning centers is important 
even at angles above the experimentally defined 9 a . 

We now discuss the value of the exponents in the field- 
(B ~ B^) and angular regime (f 5° < 9 < 65°) in which 
fully anisotropic scaling is observed. We assume that 
the anisotropic scaling is due to the different correla- 
tion length exponents in the three directions. Then, 
6c ~ £ x , Zy ~ £ and ^ ~ £ C , where £ ~ |T - T e \~ v . We 
also assume dynamic scaling for critical slowing down, 
described by the correlation time r ~ £ z , where z is 
the dynamic critical exponent. Following the procedure 
outlined in Ref. [Toj . the longitudinal resistivities are 
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FIG. 4: Phase diagram for untwinned heavy-ion irradiated 
YBCO single crystals. An effective magnetic field is used to 
account for the anisotropy of YBCO. The symbols show the 
vortex solid-to-liquid transition temperatures obtained from 
measurements of p y at constant 9 — 45° and several fields 
(■,•); and at constant poH — 2 T and different 9 (o). At 
fields far from = 2 T, isotropic glass exponents are found 
(■) and the data can be fitted to Eq. ((2} (solid line). For 
fields comparable to B$ the transition lies above this line. At 
9 ~ 15° the transition bends towards the isotropic glass line; 
nevertheless, it remains above it up to 9 ~ 65°. Inset: p y (9) 
at 1 T and 90.2 K. At 9 = the field is aligned with the 
columns. The accommodation angle 9 a is usually defined by 
the maximum of p(9) (dashed lines). 



Pi = f ~ (T - T c ) s *, where s x = v{z - 1 + X - 0, 
s y = v{z + 1 — x ~ C) an d s z = v{z — 1 — x + 0- Thus, 
the observation of three different Sj's suggests that the 
correlation length diverges with different exponents in 
all directions. An additional constraint is needed for ex- 
tracting x an d C horn the experiments. Given the simi- 
larity between the Bose glass and the transverse Meissner 
transitions flpj ] on the one hand and the transition ob- 
served here on the other, one may argue that the relation 
£ — x~ 1 = 0, which follows if the compressibility remains 
finite at the transition, should remain valid. If this is 
so, the experimentally observed Si yield v = 1.0 ± 0.4, 
X = 1.6 ± 0.7, C = 2.6 ± 0.7 and z = 6.9 ± 2.2. This can 
be compared with the results on the transverse Meiss- 
ner transition (y ks 0.70, X — 2 and £ = 3), indicating 
that the transition observed here may be of somewhat 
different nature. 

In summary, for certain angles between the field and 
the columns and at vortex densities close to the defect 
density, the resistivity of heavy-ion irradiated untwinned 
YBCO vanishes as a power law of reduced temperature 
with different exponents in all three directions. This 
shows that the scaling properties of the phase transi- 
tion are fully anisotropic in three dimensions. Thus, our 
findings provide experimental evidence for a new kind of 
universal scaling behavior at a second-order phase tran- 
sition. It would be of interest to search for anisotropic 



scaling behavior in other physical systems with a similar 
symmetry breaking as the one considered here. 
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